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S ■ Abstract 

A random multiplicative process with additive noise is described by 
a Langevin equation. We show that the fluctuation-dissipation relation 
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q ■ is satisfied in the Langevin model, if the noise strength is not so strong. 
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1 Introduction and a Langevin equation 

The Einstein relation and the fluctuation-dissipation theorem are impor- 
tant relations in the nonequilibrium statistical physics. |fj, |2] There are 
many attemps to extend those relations in systems far from equilibrium ft, 
£|. Langevin type equations have been recently studied to describe some 
■ nonequilibrium systems such as ratchet models |(| |7| and random mul- 

tiplicative processes || 10 1. The random multiplicative processes are 



C3 



O 
O 



X 



used to describe large fluctuations in chaotic dynamical systems and eco- 
nomic activity. The Langevin equation with multiplicative noise generates 
a power-law distribution in contrast to the Gaussian distribution in thermal 
equilibrium. We study the fluctuation-dissipation relation in this Langevin 
equation. 

We study a Langevin equation for a random multiplicative process with 
additive noise: 

^ = (-i + m)y(t) + v(t), (i) 
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where y(t) is a stochastic variable, 7 is a damping constant, £(t) denotes 
multiplicative noise and rj(t) addtive noise. The noises are assumed to be 
Gaussian-white noises satisfying 

mm) = 2D5{t - t% (v(tW)) = 2T8(t - t'). 

If D is zero, this process is equivalent to the Ornstein-Ulenbeck process with 
temperature T a where ksT a = T/7. The fluctuation-dissipation theorem 
and the Einstein relation are explicitly solved in this Ornstein-Ulenbeck pro- 
cess. We will show a generalized fluctuation-dissipation relation for Eq. (1) 
with nonzero D. 



2 Time correlation and the fluctuation-dissipation 
relation 

The Fokker-Planck equation for the Langevin equation (1) is written as 
9P d . n . n d ( d \ d 2 P 

- = -- { - lvP)+D - v (- vP y T — (2) 

where P(y, t) is the probability density function for y. The stationary solu- 
tion Po(y) of Eq. (2) satisfies 

8P __( S f + D)y 

dy ~ T + Dy^°- [6) 
The stationary distribution has a form 

= (T + J Dy2)(7+^)/(2D)' ( 4 ) 

where Co is a normalization constant. The normalization constant is given 
by 

_ T e / 2 D l / 2 T{l/2 + e/2) 
c ° ~ r(l/2)r(e/2) ' 

where e = j/D, and T(x) is the T-function. The variance of y can be 
calculated as 

/ 2\ r° 2 nMJ T r(3/2)T(e/2 — 1) T 

for D < 7/2. The variance diverges for D > 7/2, since Po(y) in Eq. (4) 
decays more slowly than |y~ 3 | for y — ► ±00. We have performed a numerical 
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simulation of the Langevin model (1) using the Heun method with timestep 
At = 10~ 3 or 10~ 4 to check the theoretical results. Figure 1 displays the 
numerically obtained probability distribution function (solid curve) and the 
distribution Po(y) in Eq. (4) (dashed curve) for 7 = 1, T = 1/2, D = 3/4 
in the logarithmic scale. The tail of the probability distribution obeys the 
power law with exponent 1 + e = 7/3. 
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Figure 1: Stationary probability distribution Po(y) for Eq. (1) at 7 = 1,T = 
1/2 and D = 3/4. The histogram is the numerically obtained distribution 
and the dashed line is the distribution function (4). 



The Fokker-Planck equation and the stationary distribution are related 
to the Tsallis statistics [Tl|] . If E(y) = 1/2(7 + D)y 2 is interpreted as 
generalized energy and q is assumed to be 1 + 2.0/(7 + D), the stationary 
distribution is rewritten as Pq(E) oc 1/{1 + (q — 1)E/T} q ~ 1 , which has 
a form of the equilibrium distribution in the Tsallis statistics. If O = 0, q 
becomes 1 and the Boltzmann-Gibbs statistics is recovered. The generalized 
free energy F q is defined in the Tsallis statistics as 



F q = U q - TS q 



1- f (X LP g (y)dy 
E{y)Pi{y)dy - T J ~°° , W dy, 
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where U q is the generalized internal energy and S q is the Tsallis entropy. 
The generalized free energy decreases monotonically in the time evolution 
of the Fokker-Planck equation, since 



dF^ 
dt 



OO 1 Pip 

( 1 + D)y 2 qPi- 1 —dy + 

-00 ^ Ot q 



BP 

qP^-r-dy, 
ot 
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/oo f)p 
P*-2{(T + Dy 2 ) — + (j + D)yP} 2 dy<0. (5) 
-oo oy 



This is a kind of the H-theorem in our stochastic process [12, 13]. The 
generalized free energy takes a minimum at the stationary distribution which 
satisfies Eq. (3). 

The solution of the linear equation (1) is explicitly given by 

y (t) = y(0)e-^ +s ® + f e^'-^-SCH^rKt'X (6) 
Jo 

where S(t) = f„ £(t")dt" . The time correlation c(t) = (y(0)y(t)) is given by 

C(t) = (y(0)y(t)) = (y(Of)(e-~^), (7) 

since there is no correlation between y(0) and rj(t) and the average value of 
rj(t) is zero. The probability distribution of is the Gaussian distribution 
of variance 2Dt. The time correlation is therefore given by 

C(t) = T e -(T- g )*. (8) 
w 7 - 2D v 7 

The time correlation decays exponentially and the decay constant is 1/(7 — 
D). The decay constant diverges at D = 7. For 7/2 < D < 7, the variance 
of y diverges but the decay constant does not diverge. Figure 2 compares the 
numerically calculated correlation function c(t) = (y(0)y(t)) /(y(0) 2 ) (solid 
curve) with the exponential decay e~^~ D ^ (dashed curve) at 7 = 1, T = 
1/2, D = 1/8. 

We have also investigated the response function for an external periodic 
force to Eq. (1). The model equation is given by 

jjf = (-7 + t(t))y(t) + / exp(iwf) + (9) 

where / is the strength and u is the frequency of the external force. The 
probability distribution function Pf for the static force with uj = obeys 
the Fokker-Planck equation including the external force term as 

dP t a r , , a / a \ a 2 p f , . 

The stationary solution to this Fokker-Planck equation is given by 

p _ expl/tan-^VWr^/v^P} n . 

/0 /0 (T + L>y 2 )(7+-D)/(2^) ' { > 
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Figure 2: Time correlation function c(t) for Eq. (1) at 7 = 1,T = 1/2 and 
D = 1/8. The solid curve denotes a numerically obtained time correlation 
and the dashed curve denotes e~( 7 ~ D )*. 



where cjq is a normalization constant. The average value (y) for sufficiently 
small / can be calculated as 



/oo 
dy 
-00 



y tan' 1 (^D/Ty) 



TD{T + Dy2)( 7 +D)/(2D) 

= fcoT-'PD- 1 ' 2 -?— r / \coa9y- 1 d0 = -f—, (12) 
J — D Jo J — D 

since cjo ~ Co for sufficiently small /. The static response function is simply 
written as Xst = (7 — D)~ l . 

The general response function for the periodic force /exp(za;i) is ex- 
pressed as 

(y(t)) = (x'(u)+ix"W)fe i " t , 

where x{ u ) = x'( aJ )+^x"( w ) is called the complex admittance. The Langevin 
equation is a linear equation and therefore the linear response relation is 
naturally expected. The fluctuation-dissipation relation can be generally 
expressed as 

G( W ) = ^^(y% (13) 

XstU 

where G(ui) is the Fourier transform of the time correlation function C(t). 
The standard fluctuation-dissipation theorem is obtained, if the thermal 
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equilibrium condition (y 2 ) = ksT a Xst is further assumed. The time cor- 
relation function and the static response function Xst have been already 
calculated and then the imaginary part of the response function is expected 
as 



-oj 



*» = ~ n, x ,„, (14) 



(j-D) 2 + uj 
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Figure 3: Imaginary part °f the complex admittance for 7 = 1, T = 

1/2 and D = 1/8. The numerical results are shown by points and the dashed 
curve denotes the function (14). 



On the other hand, the solution to the linear equation (9) is explicitly 
given by 

y(t) = y(0)e^ t+H M + f e 7(t'-*)-s(t')+3(*){ /e ^' + v (t')}dt'. (15) 
Jo 

The average time evolution of y(t) after a long time is given by 

ly(t))~ / ^'-t)-^')+m )f ^t> dt , = / e ( 1 -D){t'-t) f ^t' dtl !± 

Jo Jo 7 - D + iu> 

(16) 

where the relation (exp{H(t) — S(t')}) = exp{D(t — t')} is used, since the 
stochastic variable E(i) — S(i') obeys the Gaussian distribution with variance 
2D(t — t'). The complex admittance is calculated as 

'< \ l"g 
X H = ( 7 -D) 2 +^ 
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Equation (17) is equivalent to Eq. (14), that is, it is shown that the fluctuation- 
dissipation relation is satisfied in our model equation. 

We have numerically checked the fluctuation-dissipation relation. We 
have calculated a Langevin equation with a sinusoidal force: 

^ = (-J + mMt) + fcos(iut). (18) 

We have neglected the additive noise to simplify the numerical simulation. 
The imaginary part of the response function is numerically estimated from 

Trr 



—2 f lm 

y(t)sm(ut)dt, 

JJ-m JO 



where T m is a time interval to take a long-time average. The time interval 
T m = 20000 is used in our numerical simulation. Figure 3 displays the 
results for 7 = 1 D = 1/8 and / = 0.1. The points denote the numerically 
estimated values and the solid curve is the function (14). We can see that 
the fluctuation-dissipation relation is satisfied in this random multiplicative 
process. 



3 Browinian motion and the Levy flight 

The Ornstein-Ulenbeck process is the simplest model for the Brownian mo- 
tion. For the Brownian motion with velocity y(i), the position x(t) of the 
particle obeys 

|=»M. d») 

It is a problem to find the diffusion coefficient of the Brownian particle. The 
time correlation decays exponentially for D < 7. Then, the variance of the 
displacement Ax(t) = x(t) — x(0) is calculated as 

(Ax(t) 2 ) = f f\y(t')y(t")dt'dt" = 2(y(0) 2 ) f f~* c(t")dt"dt' 
Jo Jo Jo Jo 

(20) 



(7-D)(-,-2D)' 



for D < 7/2. The diffusion constant of the Brownian particle is written as 
Tj (7 — D)(j — 2D). This is a generalization of the Einstein relation between 
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Figure 4: Mean square (Ax(t) 2 ) of the displacement as a function of t for 
Eq. (1) and Eq. (19) at 7 = 1,T = 1/2 and D = 1/8. The numerical results 
are shown by the solid curve and the theoretical estimate (20) is shown by 
the dashed line. 



the diffusion constant and the mobility of the Brownian particle. Figure 4 
displays (Ax(t) 2 ) as a function of t (solid curve) and the theoretical line (20) 
for 7 = 1, T = 1/2 and D = 1/8. 

When D > 7/2, the variance of y diverges and the central limit the- 
orem cannot be satisfied, and therefore the Gaussian distribution is not 
expected for the displacement Ax. For D < 7, the time correlation de- 
cays exponentially. The displacement Ax(t) may be a summation of the 
small displacement in a time interval r of the order of the decay constant. 
The small displacements in the time interval r are expected to be indepen- 
dent stochastic variables, since the time correlation decays during the time 
interval. The small displacement is estimated as Ax(t) ~ y(t)r and the 
stochastic variable y obeys the stationary distribution with the power-law 
tail. The distribution of summation of n independent stochastic variables, 
each of which obeys the distribution with a power-law tail, approaches the 



"stable distribution" with the same power-law tail. [14] Our random walk is 
therefore considered to be a kind of Levy flight in one dimension. 

For the "stable distribution" with the power-law tail of exponent (3 where 
[3 > 3, the summation of n independent variables is expected to satisfy 

n 

«5>n~n^, (21) 
i=l 
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Figure 5: Time evolution of x(t) for Eq. (1) and Eq. (19) at 7 = 1,T = 1/2 
and D = 3/4. The large fluctuation is characteristic of the Levy flight in 
one dimension. 



where yi denotes the zth independent stochastic variable and a = j3— 1. [|l4|] 
In our model, the total time interval t is proportional to n and the following 
relation is expected. 

(Ax(t) a ) ~ t a l a . (22) 

The average (Ax(t) 2 ) will diverge, since the variance diverges, however, the 
average (Ax(t) a ) for a < 2 may be converged. We have calculated the 
average (|Ax(i)|) for 7 = 1, T = 1/2 and D = 3/4. Figure 5 displays the 
time evolution of x(t). The occasional large random jumps are characteristic 
of the Levy flight. Figure 6 displays the average {\Ax(t)\) as a function 
of t in a logarithmic scale. The average displacement obeys a power-law 
approximately. The exponent (3 of the power-law tail for the parameter 
values is expected to be 7/3 and therefore a = 4/3. The power-law evolution 
of (| Ax(t)\) is expressed as t s and the exponent S is expected to be 3/4, which 
is different from the value 1/2 for the normal diffusion. The slope of the 
numerically obtained line is approximately 0.76, which is close to 3/4 and 
definitely different from 1/2. 

4 Summary 

We have shown that the Langevin model with multiplicative noise exhibits 
the probability distribution with the power-law tail, however, the fluctuation- 
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Figure 6: Double logarithmic plot of (|Ax(i)|) vs. t for Eq. (1) and Eq. (19) 
at 7 = 1,T= 1/2 and D = 3/4. The dashed line denotes a line with slope 
3/4. 



dissipation relation is satisfied when the multiplicative noise is weak. When 
the multiplicative noise is sufficiently strong, the variance diverges and the 
Levy flight appears for the corresponding Brownian motion, whose velocity 
obeys the Langevin equation. The Langevin equation with multiplicative 
noise may be a simple and instructive model in the nonequilibrium statisti- 
cal mechanics. 
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